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Q^l decays of the seesaw Majorana neutrinos proposed by Fukugita and Yanagida. How- 

^ ■ ever, in order to avoid the wash-out of the produced lepton number we propose the 

• • . production of the Majorana neutrinos to occur non-thermally and sufficiently late. 

I After this time, in consequence, the B — L (baryon minus lepton) quantum num- 

^ ■ ber becomes a good "accidental symmetry" protecting the asymmetry produced. 

^ . This non-thermal leptogenesis at late time is realized by a boson decaying into the 

Majorana neutrinos with a long lifetime. Suggestively this boson could correspond 
to a scalar field which causes the cosmic inflation, the inflaton, and thus its decay 
means really the reheating of the Universe. We find that this mechanism works well 
even if the lightest Majorana neutrinos are not produced sufficiently or not present, 
and the decays of the heavier seesaw Majorana neutrinos can be responsible to the 
baryon asymmetry in the present Universe, as we illustrate by the example of the 
family replicated gauge group model. 
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1 Introduction 



Matter-antimatter asymmetry in the present Universe is one of the biggest puzzles in 
particle physics as well as in cosmology. This baryon asymmetry is usually expressed by 
the ratio of baryon (minus anti-baryon) number density to the entropy density sq in 
the present Universe as 

— = (3.7-8.8) X 10"" . (1) 

So 

It is only one, but mysterious, number in nature which we would like to understand. 
If our Universe experienced the inflationary stage in the beginnings of the history, the 
primordial baryon asymmetry would be diluted away and be essentially zero. The observed 
asymmetry in Eq. (|l]) should, therefore, be generated after the inflation. Such a generation 
mechanism is called baryogenesis and various scenarios have been proposed so far. 

Evidence of the neutrino oscillations gives an important clue for baryogenesis, since 
Fukugita and Yanagida proposed that lepton-number violation in nature might ac- 
count for the present baryon asymmetry. The experimental data suggest tiny but non- 
zero masses for the neutrinos. Introducing the right-handed Majorana neutrinos having 
heavy masses is the natural set up to explain such neutrino masses through the seesaw 
mechanism If this is the case, the lepton number violation which is crucial for lep- 
togenesis is naturally explained. In fact, non-equilibrium decays of Majorana neutrinos 
can produce a lepton number in the early Universe, which is partially converted into a 
baryon number through the electroweak sphaleron processes Therefore, the Fukugita- 
Yanagida mechanism p|, called leptogenesis, is probably the most attractive possibility 
to generate dynamically the observed baryon asymmetry in the present Universe. 

There is a variety of scenarios for leptogenesis in the literature @, 1^, H, 0. Here 
we restrict ourselves to leptogenesis via decays of the seesaw Majorana neutrinos {Ni: 
i = 1,2,3) having the hierarchical masses M3 ^ M2 ^ Mi. It is usually considered 
that the decays of the lightest Majorana neutrinos Ni are responsible for the present 
baryon asymmetry, although the decays of N2 and also generate a lepton asymmetry. 
This is because the A^i, having the lightest mass, can remain in thermal equilibrium 
after the decays of N2 and N3, and may induce additional rapid processes changing the 
lepton number (other than the sphaleron process). These processes wash out the lepton 
asymmetry from N2 and "too much" to make N2 and work as the producers of 



baryon number [jl0|, and hence the resultant asymmetry only comes from the decays of 
iVi. 

In this paper, however, we point out that the decays of the heavier Maiorana neutrinos 
N2 and A^3 can be a dominant source of the present baryon asymmetry if they are pro- 
duced non-thermally and also at very late time. The N2 and N3 produced non-thermally 
decay immediately after the production and generate the lepton asymmetry through the 
Fukugita- Yanagida mechanism. Furthermore, this production time is so late that the 



'^See Refs. M for early application of this mechanism. 
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cosmic temperature is low enough to prevent thermalization of the B — L (baryon minus 
lepton) asymmetry due to the hghtest Majorana neutrinos A^i as well as the wash-out 
processes caused by the existence of A^^i. As a result, after that time the B — L quantum 
number becomes a good "accidental symmetry", which ensures the sphaleron conversion 
from the lepton asymmetry from N2 and N-^ into the baryon asymmetry. This scheme 
is easily realized by a scalar field which dominates the energy of the Universe when it 
decays into the Majorana neutrinos with a long lifetime. Suggestively this boson could 
be described by a scalar field which causes the cosmic inflation, the inflaton, and thus its 
decay means really the reheating of the Universe [0. 

The leptogenesis picture by the /ieawer Majorana neutrinos produced non-thermally in 
inflaton decays works well even if the leptogenesis by the A^i decays is ineffective and also 
even if the iVi-neutrinos are absent in the early Universe. This means that our proposed 
scheme for the leptogenesis can be applied to a wider class of models. For example, we 
illustrate it in the family replicated gauge group model [0, 112 1 in which all the lepton and 
quark masses and mixing angles are well fitted order of magnitudewise. The proposed 
leptogenesis is well suited for this model, since the leptogenesis by the decays of Ni is not 
sufficient to explain the phenomenological baryon asymmetry in this model. It is found 
that the proposed leptogenesis works well even in this model and also, interestingly, we 
can determine values of the reheating temperature as well as the inflaton mass from the 
present baryon asymmetry using this family replicated gauge group model. 

The organization of this paper is as follows: In Section |^ we briefly review the leptoge- 
nesis by using the inflaton decays and we discuss the possibility of the baryon asymmetry 
generated by the decays of the heavier Majorana neutrinos. In Section ^ we review the 
family replicated gauge group model in which we illustrate our idea. In Section ^ we 
then derive the B — L asymmetry obtainable by the proposed ideas. Some problems of 
naturalness of the rather isolated inflaton or just scalar needed are discussed in Section ^ 
by suggesting some possible - speculative - solutions. Our conclusions go into Section IH. 



2 Leptogenesis in the inflaton decays 

Let us start by explaining briefly the leptogenesis in inflaton decays scenario (see Ref. 0). 
After de Sitter expansion of inflation ends, the inflaton decays, when the Hubble parameter 
of the Universe, H, becomes comparable to the decay width of the inflaton T^. The 
vacuum energy of the inflaton 0-field is completely released into decay products, and the 
Universe is reheated through their thermal scattering. The temperature at this time, the 
reheating temperature, Tr, is given using F,^ by 

Tr = 0.55s/mS~4> , (2) 

where = 2.4 x 10^^ GeV is the reduced Planck mass. 

We consider the seesaw Majorana neutrinos (i = 1, 2, 3) produced by the decays of 
the inflaton (p. If Ni are produced ^^non-thermally" , the ratio between the number density 



2 



of produced Ni and the entropy density s is estimated to be 

^=Br.|^, (3) 
s 2 

where Brj denotes the branching ratio of the decay channel (p ^ N^. When the decay 
rate of is much larger than T^, the decays immediately after being produced by the 
inflaton decays. (As we shall discuss later, we are interested in low reheating temperatures 
which ensure this condition.) 

Decays of Majorana neutrinos Ni break the lepton-number conservation and have the 
CP violation. There are two classes of decay channels: 

iVi ^ 0WS + ^ , Ni^(Pl, + e , (4) 

where (f)^s and £ denote the Weinberg-Salam Higgs and lepton doublets in the Standard 
Model, respectively. The lepton asymmetry generated by decays can be expressed by 



r(iv,^0^3 + £)-r(iv, 



_ -L yiVj ^ S^WS -T <■) — 1- <^1H yws "T '^J 



Just after the reheating completes, we obtain the lepton asymmetry induced by decays 
as 

— = 2^ei Bii . (6) 

We will assume for a while an accidental B — L conservation (in the sense that it is 
not imposed upon the model, but comes out). If the lepton asymmetry in (^) is produced 
well before the electroweak phase transition of the thermal history, i.e., ^ 100 GeV, 
the B and L conversion by the sphaleron process is active, and it brings a part of this 



lepton asymmetry into the baryon asymmetry as ||T3 

_ _28nL 

s ~ 79 s ' ^ ' 

Here we have assumed that there is only one weak Higgs doublet. Finally, we obtain the 
following expression of the produced baryon asymmetry 

= <.B..I« . (8) 

s 79 ^ 

Now we are at the point of justifying the assumptions which we have made in the 
above discussion. First, the production of Ni in the inflaton decays is available only when 

> 2Mi . (9) 
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Second, the estimation of the lepton asymmetry in Eq. is obtained under the require- 
ment that the iVj are produced non-thermally by the inflaton decays, which leads to the 
following condition on the reheating temperature 

Tr <KiM, , (10) 

where Ki are constants of order one defined by the decoupling temperatures of Ni to be 
rpdec _ Therefore, considering the inflation model satisfying the conditions (|^) and 

(p!0|), the decays of the Majorana neutrinos Ni generate the lepton asymmetry which is 
given in Eq. just after the reheating. 

However, the lepton asymmetry might be washed out, after it is produced, by the 
lepton-number violating processes. The most dangerous ones are the processes mediated 
by A^^i, since Ni is the lightest Majorana neutrino, so that it survives and still can be 
produced in the thermal bath after N2 and N3 have disappeared. If those processes are 
well in thermal equilibrium the produced lepton asymmetry is washed-out strongly |T0| . 
To avoid these wash-out processes, we have to invoke that the production of the Majorana 
neutrinos {i.e. lepton asymmetry) takes place at a sufficiently late time so that the wash- 
out processes have already decoupled and been ineffective. Thus, we have to consider 
sufficiently low reheating temperatures, 

Tr<k,Mi . (11) 

With such low reheating temperatures, the lighter Majorana neutrino(s) are completely 
decoupled from the thermal bath of the Universe, and the B — L becomes a good "acci- 
dental symmetry" for T < Tr. The conditions Eqs. (|^) and (0) justify our estimation of 
the baryon asymmetry in Eq. (P), by explaining why there is no wash-out effect. 

It should be noted that Eq. ( [TT| ) ensures the non-thermal production of the heavier 
Majorana neutrinos N2 and if > 2M3. This means that the decays of A''2 and can 
be dominant sources of the present baryon asymmetry (if 62 Br2, 63 Br3 ^ ei Bri). Further, 
this mechanism works even if there is no A^i, i.e., if Bri = 0. This feature is completely 
different from the conventional thermal leptogenesis |^ where the lepton asymmetry is 
generated by the decays of N2 (N^) at the temperature of T ~ M2(M3) 3> Mi, and 
hence the produced lepton asymmetry may be easily washed out and the resultant baryon 
asymmetry comes from the decays of the lightest Majorana neutrinos A^^i. 

These observations lead to that our proposed scheme for the leptogenesis can be applied 
to a wider class of models of the sort explaining fermion masses and mixings. In the next to 
next section, we will illustrate it in a specific example, a model presented in the following 
section. 

3 Family replicated gauge group model 

In this section we review briefly the family replicated gauge group model [0, |12|. This 
model is based on a large gauge group being the Cartesian product of family specific 
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gauge groups, namely, 

J^,.SSMG,xU{l),_,,) , (12) 

where SMd denotes SU{?>)i x SU{2)i x f/(l)i (Standard Model gauge group), Ug-L^i is 
B—L gauge group, and i denotes the generation. This group in Eq. (|12]) is represented only 
by representations that do not mix the different irreducible representation of the Standard 
Model and it breaks spontaneously down to the group SU{?>) x SU{2) x U{1) x U{1)b_l 
at the scale about 1 to 2 orders of magnitude under the Planck scale. The breaking to the 
Standard Model gauge group is supposed to occur by means of five Higgs fields which we 
have invented and denoted by the symbols W , T, p, oj and 0ss- Their quantum numbers 
are given in Table 1. Finally the breaking of SU{2) x U{1) of the Standard Model is 
caused by a Weinberg-Salam Higgs field, 0ws (also its quantum numbers are found in 
Table 1). 

We summarize here the vacuum expectation values (VEVs) of the six Higgs fields 
which the model contains: 



1) The smallest VEV Higgs field is the Weinberg-Salam Higgs field, 0ws; with the VEV 
at the weak scale, 246 GeV. 

2) The next smallest VEV Higgs field, called 0ss; is also alone in its class and breaks the 
"diagonal B — L" gauge group U{1)b_l, common to all the families. This symmetry 
is supposed to be broken (Higgsed) at the seesaw scale as needed for the neutrino 
oscillation scale in seesaw mechanism. This VEV turns out of the order of 10^^ GeV. 

3) The next four Higgs fields are called W, T, p and u and have VEVs of the order of a 
factor 4 to 50 under the Planck unit. It means that if intermediate propagators have 
scales given by the Planck scale, as we assume, they will give rise to suppression 
factors of the order 1/4 to 1/50 each time they are needed to cause a transition IT^ . 



The quantum numbers of the 45 well-known Weyl particles and the three additional 
Majorana neutrinos to be used as seesaw neutrinos are strongly restricted by the re- 
quirement that all anomalies (gauge and mixed ones) shall vanish even without using 
Green-Schwarz anomaly cancellation mechanism |15|. Thus the family replicated gauge 
group model is an anomaly free model. All the U{1) quantum charges in this model can 
be found in Table 1. 

In the family replicated gauge group model there exist many bosons and fermions at 
the fundamental scale (the Planck scale). However, in order to have no mass protection, 
the left-handed fermions and its right-handed partner must combine to a Dirac or vector- 
like-coupled fermion, in the sense that they are described as Dirac particles from the 
picture of the Weyl particles: they are combinations of left-handed and right-handed 
states with the same (gauge) quantum numbers. The left-over Weyl particles (in other 
word those without chiral partners) in this model - i.e. those that could get small masses 
relative to the Planck scale - are specified in more detail and are actually assumed to 
form a system of three proto-generations, each consisting of the 16 Weyl particles of a 
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usual Standard Model generation plus one seesaw particle. In this way we can label 
these particle as proto-left-handed or proto-right-handed w-quark, d-quark, electron etc. 
To get the quantum numbers in the model under the gauge group for a given fermion 
proto-irreducible representation, we proceed in the following way: We note the generation 
number of the particle for which we want the quantum numbers and we look up, in the 
Standard Model, what the quantum numbers are for the irreducible representation in 
question and what the B — L quantum number is. For instance, if we want to find the 
quantum numbers of the proto-right-handed strange quark, we note that the quantum 
numbers of the right-handed strange quark in the Standard Model are: weak hypercharge 
y/2 = —1/3, singlet under SU{2), and triplet under SU{3), while B — L is equal to the 
baryon number = 1/3. Moreover, ignoring mixing angles (as we do for proto quarks), the 
generation is denoted as number i = 1. That it belongs to generation i = 1 means that 
all the quantum numbers for SMGj, j = 2, 3 are trivial. Also the baryon number minus 
lepton number for the proto-generation number two and three are zero: only the quantum 
numbers associated with proto-generation one are non-trivial. Thus, in our model, the 
quantum numbers of the proto-right-handed down quark are 1/2/2 = —1/3, singlet under 
5f/(2)i, triplet under 5f/(3)i and {B - L)i = 1/3. 

For each proto-generation the following charge quantization rule applies 

1 + 1 + 1 = (modi) , (13) 

where and di are the triality and duality for the i'th proto-generation gauge groups 
SU{3)i and SU{2)i respectively, of course a consequence of the corresponding rule in the 



Standard Model (See Refs. [[11], [12 



Combining Eq. (p!3| ) with the principle of taking the smallest possible representation of 
the groups SU{?>)i and SU{2)i, is sufficient to specify the non abelian quantum numbers 
in terms of the six Abelian quantum numbers yi/2 and {B — L)i. Using this rule we easily 
specify the fermion representations as in Table 1. 



3.1 Mass matrices 

We can easily evaluate using the system of quantum numbers given in Table 1 the numbers 
of Higgs field VEVs of the different types needed to construct the transition between the 
left- and right-handed Weyl fields corresponding to a given element in a mass matrix. 

The main point is that the order of magnitudes of the mass matrix elements are 
determined by the number and order of magnitudes of the Higgs field VEV factors needed 
to accomplish a given quantum number transition. We have namely assumed that all the 
Yukawa couplings, and also other couplings, as well as the masses of particles that are 
not mass-protected are of order unity in some fundamental scale (which we usually take 
to be the Planck scale). Thereby of course the whole calculation can only give order of 
magnitudewise results. 

We shall now write down the mass matrices - five Dirac type and one Majorana 
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Table 1: All U{1) quantum charges in the family replicated model. The symbols for the 
fermions shall be considered to mean "proto" -particles. Non-abelian representations are 
given by a rule from the abelian ones (see Eq. (|13])). 
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neutrino mass matrix - which are necessary to fit the fermion masses and mixing angles 
and also to discuss the mechanism of the baryogenesis: 



the up-type quarks: 
the down-type quarks: 



/ (u^fW^T^ up^W^T^ up^{W^)^T 



V2 



t^4 



s) 



V2 



\u^pW\T^Y W\T^Y WT 



(14) 



(15) 



the charged leptons: 

(0ws) 



cu6(pt)3i^(Tt)2 WiT^y W\T^)^ I (16) 

\co%p^fiW^yT^ (W^yr^ WT 



V2 

the Dirac neutrinos: 

(i<P^sV)( (^t)3^3^tT2 (^t)3^3^2(yt)7 

/ (^t)3^tT2 w^T^ W\T^y \ (17) 
and the Majorana (right-handed) neutrinos: 

/ (pt)6r6 (pt)32^6 (pt)3^3(2^t) 

Mr ~ (0ss) (P^)'r6 M/3(Tt)3 

V(p^)W3(Tt)3 W^3(rt)3 H^6(Tt)12 
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We know neither the Yukawa couphngs nor the precise masses of the fundamental 
fermions, but it is one of our basic assumptions of naturalness of the model that these 
couplings are of order unity and random complex numbers at the Planck scale. In the 
numerical evaluation of the consequences of the model we explicitly take into account these 
uncertain factors of order unity by providing each matrix element with an explicit random 
number factor Xij with a distribution, so that its logarithm average (logAjj) ~ and its 
spread is 64%. At the end we then average the obtained results over these random number 
distributions. Note that the random complex order of unity factors which are supposed 
to multiply all the mass matrix elements are not written explicitly in Eqs. (14)- (18) but 
are understood to be there anyway. 

The philosophy of the model is that these mass matrices correspond to effective Yukawa 
couplings to be identified with running Yukawa couplings at the fundamental /Planck scale 
for the Higgs field 0ws in the case of the first three mass matrices and for (pss in the case 
of the right-handed neutrino mass matrix. Therefore these effective Yukawa couplings 
have in principle to be run down by the beta-functions to the scale of observation, see 
section |3.2| . It is also important that we include the "running" of the irrelevant operator 
of dimension 5 giving the neutrino oscillation masses. 

The right-handed neutrino couplings - or equivalently mass matrix elements - are used 
to produce an effective mass matrix for the left-handed neutrinos which we after the five 
dimensional running down mentioned take as the ones observed in neutrino oscillations 0, 

Meff~M,^M^^(M,^)^. (19) 



3.2 Renormalization group equations 

It should be kept in mind that the effective Yukawa couplings for the Weinberg- 
Salam Higgs field, which are given by the Higgs field factors in the above mass matrices 
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multiplied by the understood order of unity factors (taken as random numbers), are the 
running (effective) Yukawa couplings at a scale very close to the Planck scale. Thus, 
in our calculations, we had to use the renormalization group /3-functions to run these 
couplings down to the experimentally observable scale, i.e. fi = 1 GeV where fi is the 
renormalization point. This is because we took the charged fermion masses to be compared 
to "measurements" at the conventional scale of 1 GeV. In other words, what we take as 
input quark masses are the current algebra masses, corresponding to running masses at 
1 GeV. We used though for the top quark the pole mass instead [0: 

M.=m,(M)(l + l^^ , (20) 

where we set M = 180 GeV as an input, for simplicity. 

From the Planck scale down to the seesaw scale or rather from where our gauge group 
break down to SMG x U{1)b^l we use the one-loop renormalization group running of 
the Yukawa coupling constant matrices, Yj^, Yq, Yg, and Yj, (being proportional to the 
mass matrices M^, M^, M^, and M^, respectively), and the gauge couplings |]IB|, |TT 



rft 10^^ ' ^ ^ 

16vr^^ = -^-^gl , (22) 

16n^^ = -79l , (23) 

le^r^^ = I {YuiYuy - Y'.(r.)t) + { - [^gl + \gl + 8gi) } , (24) 

16vr^^ = I {YniY.y - Y,{Y,y) + { - [Igf + \gl + Sgj)] Y, , (25) 

^^"'^ = I (^^^^^)' - ^^(^^)') + { - (i^i + l^l)} y. , (26) 

ISvr^^ = I iYXY.)^ - r.(F.)t) Y^ + \Ys- H-gl + ^gi]] K , (27) 



levr^^ = ((1;)^K) + ({Y^)^Y:f , (28) 
Ys = Tr(3FjF^ + 3F,ty^ + yjy^ + ytyj , (29) 

where t = In /x and /i is the renormalization point. 

However, below the seesaw scale the right-handed neutrino are no more relevant and 
the Dirac neutrino terms in the renormalization group equations should be removed, and 
also the Dirac neutrino Yukawa couplings themselves are not accessible anymore. That 
means that, from the seesaw scale down to the weak scale, the only leptonic Yukawa 
/5-functions should become as follows: 

16vr^^ = \ [Y,{Y,y) Y, + [Ys- {^/, + ^,,^)} . (30) 
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Note that the quantity, Yg, must be also changed below the seesaw scale: 

Ys = Tt{3Y^Y^ + 3Y^Y, + Y^Y,) . (31) 

In fact we stopped the running down according to formula P7| ) differently for the 
different matrix elements in the Y^, matrix corresponding to the right-handed neutrino 
mass supposed most important for the matrix element in question. 

Starting the running in an analogous way, we further should evolve the effective neu- 
trino mass matrix considered as two Higgs two fermion irrelevant - a five dimensional - 
term [1^ from the different right-handed neutrino masses to the weak scale (180 GeV) 



depending on the terms: 



IGtt^^^ = i-Sgl + 2X + 2Ys) - - [M,s (Y^Y^) + {Y^Y^f M,s) , (32) 



where Yg defined in Eq. (|3lD and in this energy range the Higgs self-coupling constant 
running equation is 

levr^^ = 12A^ - {^Ig! + 9gl) A + ^ (|^,^ + ^glgl + gl) + 4F,A - 4Hs , (33) 

with 

Hs = Tr |3 {Y^Yuf + 3 {Y^f + {y^Y,)'] . (34) 

The mass of the Standard Model Higgs boson is given by Mfj = A {(p-ws)"^ and, for 
definiteness, we took Mh = 115 GeV at weak scale. 

From 180 GeV down to 1 GeV - experimental scale (1 GeV) - we have evaluated the 
beta-functions with only the gauge coupling constants. In order to run the renormalization 
group equations, we use the following initial values: 



m 

SU{2) 
SU{3) 



gi{Mz) = 0.462 , (7i(Mpianck) = 0.614 , (35) 
g2iMz) = 0.651 , (72(Mpi,nck) = 0.504 , (36) 
gs{Mz) = 1.22 , (73(Mpianck) = 0.491 . (37) 



Note that we have ignored the influence of the B — L gauge coupling constants; however, 
this effect should not be significant, because there is from the Planck scale to the seesaw 
scale (~ 10^^ GeV) only a factor lO'^. So it should be good enough for our order magnitude 
calculations. 



3.3 CP violation in Majorana sector 



The CP violation in the Majorana neutrino decays in Eq. (^, e^, arises when the 
effects of loops are taken into account, and at the one-loop level, the CP asymmetry 
comes both from the wave function renormalization and from the vertex correction 




(MP 



(38) 
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where can be expressed through the unitary matrix diagonahzing of the right-handed 
neutrino mass matrix Vr. 



MP 



M^Vr 



vIMrMIVr = diag (M,^M2^M 



(39) 
(40) 



The functions in Eq. (|38|), f{x) and g{x), are obtained by the calculations in pertur- 
bation theory, the one-loop vertex contribution and the self-energy contribution, respec- 
tively, under the condition that the differences between Majorana neutrino masses are 
sufficiently large compared to their decay widths. Their precise forms are as follows: 



fix) 



1 + x) In 



1+x' 



X 



X 



(41) 



We have taken into account the renormalization group running effects in all sectors, 
not only the charged sectors, but also Dirac- and Majorana neutrino sectors to evaluate 
from the Planck scale down to the scale, 1 GeV. Furthermore, we have even taken into 
account the running of the dimension five operator involving two 0ws's and two left- 
handed neutrinos which give the neutrino oscillation masses. The e^'s were obtained at 
the corresponding temperature, i.e., the Majorana neutrino masses (T = Mi). 



3.4 Mass and mixing results 



The calculation using random numbers and performed numerically was used to fit the 
masses and mixing angles to the phenomenological estimates by minimizing what we call 
"goodness of fit" , 

' m fit ' 



g.o.f. = 



In 



I ,exp ^ 



(42) 



a kind of for the case that we have only order of magnitude accuracy. The results 
presented in Table 2 were obtained from the following values of the set of Higgs VEVs 
- where the Higgs field VEVs for the fields W, T, p and u causing the breaking to the 
diagonal subgroup SMG x U{1)b-l are quoted with the VEV in Planck units, while they 
are for (pss (and the not fitted (pws) given in GeV units: 



{W) = 0.157 , (T) = 0.0804 , (p) = 0.265 , 
(0ss) = 5.25 X 10^^ GeV , (0ws) = 246 GeV 



{u) = 0.244 



(43) 



The results of the best fit, with the VEVs in Eq. (0), are shown in Table 2 and the fit 
has g.o.f. = 3.51. To see a typical error, say average error, compared to the experimental 
values we should divide this value with the number of predictions (17 — 5 = 12) and 
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Table 2: Best fit to conventional experimental data. All masses are running masses at 
1 GeV except the top quark mass which is the pole mass. Note that we use the square roots 
of the neutrino data in this Table, as the fitted neutrino mass and mixing parameters, in 
our goodness of fit (g.o.f.) definition, Eq. (^2p. 





Fitted 


Experiment 




Fitted 


Experiment 


niu 


4.9 MeV 


4 MeV 




4.8 MeV 


9 MeV 


rric 


0.68 GeV 


1.4 GeV 


Ills 


325 MeV 


200 MeV 


Mt 


202 GeV 


180 GeV 


rrib 


6.0 GeV 


6.3 GeV 


rUe 


1.7 MeV 


0.5 MeV 


Vus 


0.11 


0.22 




122 MeV 


105 MeV 




0.028 


0.041 


rrir 


1.53 GeV 


1.78 GeV 


Vub 


0.0030 


0.0035 


tan^ 6q 


0.27 


0.38 


Am| 


9.0 X 10-5 eV^ 


6.9 X 10-5 eV^ 


tan^ ^atm 


0.55 


1.0 




1.8 X 10-3 eV^ 


2.5 X 10-3 eV^ 


tan ^chooz 


2.57 X 10-2 


< 2.6 X 10-2 


g.o.f. 


3.51 






Prediction 


Experiment 




Prediction 


Experiment 


Jcp 


4.4 X lO"*' 


(2 - 3.5) X 10-5 


(heel) 


3.3 X 10-3 eV 


(0.11-0.56) eV 


T (p-.7r0e+) 


> 10^3 yrs 


> 2.6 X 10^3 yrs 


Br{^.^e7) 


~ 7.0 X 10-56 


< 1.2 X 10-11 



then take the square root of it: y3. 51/12 = 0.54. This means that the 12 degrees of 
freedom have each of them a logarithmic deviation of about 55%, i.e., we have fitted all 
quantities with a typical error of a factor exp (^±^3.51/12^ ~ LOOtgisg. This agrees with 
theoretically predicted deviations [|^. 

Note that to make the best fit values of the mass squared difference and mixing 
angle parameters presented in Table 2 we used the following point as fit values from the 



combination of the latest Super-Kamiokande and SNO experimental results |20 



Ami = 6-9 X 10-5 eV^ and tan^ Oq = 0.38 . (44) 

The atmospheric neutrino parameters are the following according to the Super-Kamio- 
kande results PT| : 

Aml^^ = 2.5 X 10-^ eV^ and tan^ ^^tm = 1-0 . (45) 



With the VEVs in Eq. (^3|) we can calculate several physical observables [y - Jarl- 
skog triangle area Jcp, proton lifetime t{p Tr^e"*"), the effective electro-neutrino 
mass {\mee\) relevant for the neutrinoless double beta decay p3|, and the branching ratio 
Br(/i — s> 67) - which can be found also in Table 2. 

We also present our predicted hierarchical left-handed neutrino masses (m^) and the 
right-handed neutrino masses (Mj) as well as the CP violation in the Majorana neutrino 
decays (e^) with mass eigenstate indices {i = 1,2,3): 

mi = 1.3 X 10"^ eV , Mi = 1.4 x 10^ GeV , |ei| = 1.8 x 10-^° , 
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ms = 9.6 X 10-3 eV , M2 = 1.4 x 10^° GeV , jeal = 2.2 x 10"^ , (46) 
ma = 4.3 x 10~^ eV , M3 = 1.6 x 10^° GeV , leg] = 1.8 x 10^^ . 

Due to the philosophy that all coupling constants are of order unity and complex 
random numbers at fundamental scale, i.e., the phases are rotated randomly. So we are 
not able to predict the sign of these quantities. Therefore, we present absolute values for 
the CP-violating parameters in the Majorana sector. 

We estimate for later calculations of baryon asymmetry in this model the uncertainties 
for the calculated physical quantities. The order unity complex numbers are given by a 
Gaussian distribution with mean value zero, which leads to a fluctuation in the logarithm 
of ±64%. So the expected fluctuation in the natural exponent is 64% ■ for the e's 
(straight forward assumption of independence of 0{1) factor fluctuations in numerator 
and denominator of Eq. ( pH]) would lead to 64% ■ vTT, but crude expectation of compen- 
sating correlations may justify 64% ■ a/S) and 64% for the right-handed neutrino masses, 
respectively. Here are these quantities with the errors: 

Ml = l.Atlie X 10^ GeV , |ei| = 1.8^?| x 10"^° , 

M2 = l.Atlie X 10^° GeV , |e2| = 2.2tl's x 10"' , (47) 

M3 = 1.6^^:^6 X 10^° GeV , legl = 1.8^?! x 10"^ . 



4 Estimation oi B — L production 

We shall discuss here the B — L production in the model described above. It is found 
from the results given in Eq. ( ^61 ) that the leptogenesis induced by the lightest Majorana 
neutrinos A^^i cannot account for the present baryon asymmetry (see Eq. ([l|)), since the 
decay rate asymmetry |ei| for A^^i is too small. Therefore, we have to investigate the 
possibility of the leptogenesis by the heavier ones, N2 and N^, being then produced non- 
thermally in infiaton decays. 

With well-suited assumptions on the infiaton mass and the reheating temperature 
Tr the formula (|^) gives the baryon asymmetry produced by decays of N2 and as 

^==^^i:nb.i^ ,48) 

s 79 ^<t> 

T 

^ fl.2 X 10-^Br2 + 0.96 X IQ-^Brg) — . (49) 

As we have already mentioned that the applied model is not able to predict the sign of 
the considered quantities, e's, we will use Eq. (|48|) instead of Eq. (||) for baryogenesis 
calculation. Here and hereafter we use the central values of physical quantities given in 
Eq. (E6D in the calculation. 
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First of all, we make the rough estimate on the possible value of the baryon asymmetry 
generated by the proposed mechanism. The condition on in Eq. (P) leads to 



< 1.2 X 10-^ Bra + 0.96 x 10'^ Brs ^ . (50) 
Moreover, the reheating temperature should be low enough to avoid the wash-out pro- 



cesses by A^^i as explained in Eq. ([Tl]). By taking Ki = 1, we find the produced baryon 
asymmetry to be 

— < f 1. 2 X 10"^ Brs + 0.96 X 10"^ Bra) ^ 
s ^ ^ 2M3 

= 5.1 X 10""Br2 + 4.2 X lO^^^Brg , (51) 



which gives, in the extreme case of Br2 = Bra = 1/2, 



'''' <4.7x 10-" . (52) 



s 



Therefore, broadly speaking, in the family replicated gauge group model described in the 
previous section the observed baryon asymmetry (|I|) can be explained by invoking the 
proposed non-thermal leptogenesis. 

This non-thermal scenario for the leptogenesis requires the reheating temperature in 
the region of 

Tr < Ml ~ 10^ GeV . (53) 

It is important to mention that such low reheating temperatures can be naturally obtained 
in the considered model. We expect naturally that the infiaton is totally singlet under 
the gauge groups of the model given in Eq. ( p!2[ ) in order to prevent the large radiative 
corrections disturbing the flatness of the infiaton potential. Then, we obtain the couplings 
between infiaton and Majorana neutrinos Ni as 

C = (gauge invariant terms) 

D ^c,M,N,N, , (54) 

where we have taken the effective cut-off scale as the reduced Plank scale M^. It is crucial 
to notice that the strength of these interactions are determined up to (9(1) factors Cj by the 
charges of the Majorana neutrinos under the flavor (gauge) symmetries [Q. The infiaton 
might couple to other fermions and/or Higgs fields similar to Ni, which might disturb our 
leptogenesis by raising the reheating temperature and by lowering the branching ratios 
Brj (See the discussion in the next section.). 

It is found from Eq. (Q) that the dominant channels of the infiaton decay are (p ^ 
N2N2 and/or (p ^ N^N^ when > 2M3, and the total decay rate of the infiaton is 
given by 



2\ 3/2 

_L„„,g|,f^^^J h_I- . (55) 



1 , ,2 rMi\' 4M, 
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Figure 1: The produced baryon asymmetries ub/s are shown by the (red) thick-sohd 
hues for |c2,3| = 0.1, 1 and 10, respectively. The dashed hues represent the regions where 
the reheating temperature is higher than the hghtest Majorana mass Tr > Mi. The 
typical effect of the statistical error of the CP asymmetries €2,3 is also represented for 
rrifp = 1 X 10^^ GeV and |c2,3| = 1. The observed ub/s are denoted by the region between 
the horizontal dot-dashed lines. 



The reheating temperature is estimated from Eq. (^, by neglecting the factor (1 
and by taking |c2| = |c3| = 1, as 



4M2/m2)3/2 



/M2\^ /M3\^ , 



3.1 X 10^ GeV X 



10" GeV 



1/2 



(56) 



Interestingly, it shows that the required reheating temperatures in Eq. (^) are naturally 
obtained when m,^ ~ 2M2,3 with couplings |c2,3| of order one. Therefore, the interaction 
Lagrangian ( ^If ) ensures the late decay of the inflaton, i.e., the non-thermal leptogenesis 
at late time. 



In Figure || we show the produced baryon asymmetries (|48|) by using the exact formula 
for the inflaton decay width in Eq. (0). We neglect the region where the reheating 
temperature is higher than the lightest Majorana neutrino mass Tr > Mi. It can be seen 
that the observed baryon asymmetry is obtained when the inflaton mass is close to 2M2,3 
with |c2,3| ~ 10. The inflaton mass of ~ 2M2,3 but m,^ > 2M2,3 is crucial to have 
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a sufficiently low reheating temperature Tji <, Mi in order to avoid the wash-out of the 
produced asymmetry by A^^i, as well as to have a larger baryon asymmetry. 

One might worry that we only get agreement with the baryon asymmetry for the value 
10 of the order one constants, C2 = C3 = 10, given in Eq. (^41). However, it is basically not 
a problem at all: we have used the best fitted values for the CP-asymmetry parameters 



e2,3, but they have large statistical uncertainty as represented in Eq. (|47D . We show this 
typical error also in Figure |1|. It is found that the observed baryon asymmetry is easily 
obtained with |c2,3| ~ a few, which is very consistent with the philosophy that every 
dimensionless coupling takes an order one value at the fundamental scale. 

To summarize, we have found that the proposed non-thermal scenario for the lepto- 
genesis indeed works well in the family replicated gauge group model, and the dominant 
contributions to the present baryon asymmetry comes from the out-of-equilibrium decays 
of the heavier Majorana neutrinos A''2 and rather than the lightest one, A^i. To avoid 
the strong wash-out of the produced lepton asymmetry mediated by Ni, it is crucial that 
the infiaton decays into N2 and at the late time via the interaction terms in Eq. (0), 
which ensures naturally the reheating temperature of sufficiently low. Interesting enough, 
we have observed that the successful baryogenesis requires (or predicts) the mass of the 
infiaton of ~ 2M2,3 ~ 10^° GeV. 



5 Discussion 



5.1 Problem of decay into Weinberg-Salam Higgs particles 

For the successful baryogenesis the infiaton particle should survive until the era of the 
temperature T<,Mi, when the "accidental" conservation law of the B — L charge has 
been installed to prevent the wash-out of the produced lepton asymmetry. Furthermore, 
the branching ratio of the infiaton decays into the heavier Majorana neutrinos N2 and 
iVs should be sufficiently large to account for the observed baryon asymmetry. To realize 
this long-lived infiaton it is crucial that the total decay width of the infiaton is really 
dominated by the decays into the Majorana neutrinos N2 and A^3 via the interaction (0). 

In the family replicated gauge group model discussed above, other infiaton decay 
channels than the Majorana neutrinos do exist. In fact, the infiaton can also decay into 
pairs of quarks or charged leptons. These decay amplitudes are proportional to the masses 
of the produced fermions, and hence it is naturally understandable that the heaviest 
fermions N2 and A^3 will dominate among the final states of two fermions. Moreover, by 
assuming the philosophy that the dimensionless couplings are of order one, we naturally 
expect that the Higgs fields in the considered model obtain the masses of order of their 
VEVs. Among these the Higgs fields, the masses of W , T, p, uj, and 0ss are sufficiently 
heavy (see Eq. (|i3|)) that the infiaton of mass ~ 2M2^3 cannot kinematically decay 
into them. The infiaton also decays into pairs of gauge bosons in the Standard Model, i.e., 
photons, ly^-bosons, Z°-bosons, and gluons. These decay channels would be potentially 
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dangerous depending on the inflaton mass. However, in the inflaton mass region we are 
interested in {i.e., ~ 2M2,3 but well above threshold) the partial widths of decays into 
the Standard Model gauge bosons are comparable to those of — >■ A^2,3-^2,3- Therefore, 
they do not disturb our successful leptogenesis scenario. Other gauge fields appear in the 
considered model obtain the masses from the Higgs fields W , T, p, uj, and (pss, and the 
inflaton decays into these gauge fields are kinematically forbidden. 

However, there is one point that needs discussions: why should the inflaton not decay 
into a couple of Weinberg-Salam Higgs bosons, one Higgs and one anti-Higgs? A priori 
there could exist interaction terms in the Lagrangian density of the form 

£ = ;U |0ws|^ + |0ws|^ , (57) 

where /i is a mass parameter and A a coupling constant. These interactions allow the 
inflaton decay with a width 

r(0 0WS 0ws) = Q • (58) 

That would, however, have the adverse effect of allowing the fast decay of the inflaton 
and making the formal reheating temperature Tr in Eq. (^ higher than the mass of (p, 
m^, since we expect from our philosophy the coupling A be the order unity and also the 
inflaton and its VEV be of the same order, ~ (0). It means that the inflaton 
would decay immediately after the inflation, with exceedingly small branching rations to 
seesaw neutrinos and our picture would be spoiled. In practice, in order that the 0ws 0ws" 
channel leave the Majorana neutrino channel(s) to dominate, the required suppression of 
the effective coupling A^^ compared to unity is found from Eqs. (|55|) and (|58D to be of the 
order of 

- |c.3| ^ - 10 . |C,3| [j^^) , (59) 

where we assumed that the mass of the inflaton and its VEV are of the same order, 
rrifj,/ (0) ~ 1. Therefore, to realize the successful leptogenesis the coefficients A and /// (0) 
should be extremely small contrary to our basic philosophy. This is a drawback of the 
proposed scenario. 

However, first of all, we must say that we had already accepted this kind of finetuning 
of parameters in the model. Since we have no symmetry to control radiative corrections, 
we would expect every mass parameter in the model to be of the order of the fundamental 
cutoff scale, e.g., the Planck mass, under the assumption of dimensionless coupling being 
order unity. On the contrary, we had already introduced the mass scales which is smaller 
than the Planck scale, which are the VEVs (or masses) of the Higgs fields. The most 
severe one is the mass of the Weinberg-Salam Higgs. Although the hierarchy in various 
Yukawa couplings to explain the observed mass spectrum of fermions can be very elegantly 
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explained by the flavor (gauge) symmetry, we had not explained the hierarchy between 
the Planck scale and the mass scales of interest, but just realized it by finetuning by 
hand. In fact, as an example, the order unity coupling X^g = C(l) also induces the huge 
correction to the mass of the Weinberg-Salam Higgs particle, i.e., Srriy^s = VKs (0) ~ 
a/ Aeft ~ 10^° GeV. To have the weak scale mass ~ 10^ GeV for the Weinberg-Salam 
Higgs, the coupling would be extremely small X^g ^ 10^^^, which is more stringent than 
Eq. (|59|). However, this is not a matter at all from our standpoint, since, in addition to 
other dangerous (radiative) corrections, such a huge contribution is tremendously tuned 
so as to have the Weinberg-Salam Higgs mass of the weak scale ~ 100 GeV. In the 
considered model we had already accepted the finetuning for the mass parameters in the 
scalar sector. Therefore, the required suppression of the coefficients as given in Eq. (|59D , 
which makes the non-thermal leptogenesis works well, is a problem of our attitude toward 
the finetuning. Although we loose the simple philosophy of couplings being order unity, 
the successful leptogenesis might be obtained just by accepting not only the finetuning 
for the hierarchy problem, but also of the dimensionless couplings appearing in the scalar 
sector. Moreover, it should be remarked that infiaton scheme in general suffers from 
difficulties to allow the infiaton interactions with the Standard Model particles for the 
reheating, since the interactions should be sufficiently weak not to spoil the slow-roll 
infiation scenario. 

Although the required suppression in Eq. (|59D would be easily achieved by the finetun- 
ing, hereafter we will briefiy discuss two possible resolutions of the finetuning problems. 



5.1.1 Multiple point principle "explanation" of suppression of infiaton decay 
to Weinberg-Salam Higgses 

One possible attempt to explain the absence of the terms, e.g., X |0ws|^ 0^ term consists 
in calling upon a principle which has actually been developed in connection with the 
model discussed here as the example. This is the postulate of the multiple point principle 
(MPP) saying that the effective potential as a function of the scalars such as 0ws and 
(f) etc. should have many (as many as possible) equally deep minima. This is a principle 
that is actually true in superdsymmetric models in as far as there all the superdsymmetric 
states have to have zero and (thus) minimal energy density. 



In the article it was argued that such a principle of equally deep minima would 
have a tendency to split into separate sectors - totally decoupling ones - once there are 
only one (or a few) free parameters adjustable interactions between potentially separable 
sectors. 

Ignoring the irrelevant terms such as Eq. (^) the only interaction between an "infiation 
sector" and "our sector" - meaning the Standard Model fields plus the seesaw neutrinos 
as well as possible gauge fields or scalars associated with it - is supposed to be the term 
A |0ws|' 101' in Eq. (0). 

Now we want to argue that provided we could find solutions for the coupling constants 
with, say, two degenerate minima in both of the mentioned sectors, when they were 
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separate, we can argue for that A = is at least one possible solution for satisfying the 
MPP requirement. Let us argue for that by counting relations (~ finetunings) between 
the coupling constants needed to achieve degenerate minima in effective potentials: to 
achieve n degenerate minima in an effective potential V^s - it be a function of one or 
more scalar fields - one needs to finetune {n — 1) couplings or parameters. Let us for 
example imagine that if we had our sector alone MPP could tune in 2 — 1 = 1 relation 
between the couplings and thereby achieve 2 degenerate minima (in the Weinberg-Salam 
Higgs effective potential) and that one also by tuning 2 — 1 = 1 coupling/parameter can 
achieve two generate minima in the infiaton effective potential, if that were alone. Then 
we can see easily that by combining the two sets of coupling constants and taking further 
A = 0, we achieve 4 = 2-2 degenerate minima for the effective potential of the full model. 
This effective potential is (at least) a function of both and and it has its degenerate 
minima in all the (0, 0ws)-pairs which are obtained by combining the 0-values from the 
degenerate minima in the 0-potential alone and the 0ws-values from the generate minima 
of "our sector" alone. There are 2-2 = 4 such combinations (0, 0ws) and thus this special 
solution with A = provides 4 degenerate minima. A priori we expect that 4 degenerate 
minima requires 4 — 1 = 3 finetuning relations between the couplings/parameters. The 
proposed solution with A = uses just 3 finetunings in as far as we used one finetuning 
relations for each of the two sectors plus the finetuning A = 0. The proposed solution 
has therefore not used more finetunings to be fixed by MPP than what is expected to 
be needed anyway. Had we instead imagined that we had got arranged say 3 degenerate 
minima in the two sectors separately we would have been able to produce a solution with 
A = and 9 generate minima, that would have needed only 2 + 2 + 1 = 5 finetunings 
against the a priori expected 9 — 1 = 8 in this case. So in this case we would quite clearly 
get A = 0. We see here the possibihty for the MPP to produce decouphng (totally) of 
only loosely coupled sectors, such as the infiaton one and ours. 

The region in field space in which this kind of argument fixes A conceived of as renor- 
malization group running with the fields (0, 0ws) is of course where |0|^ |0ws|^ is large. 
So a priori one may wonder if the running of A could let it be non-zero for some other 
field values. It is, however, rather easily seen that the /3-function for A only obtain terms 
proportional to A itself. Thus once zero the coupling A will remain zero under the running. 

5.1.2 Supersymmetric Extension 

One of the most natural ways to explain the huge mass hierarchy between the fundamen- 
tal Planck scale and the electroweak scale is to introduce an additional symmetry, super- 
symmetry (SUSY). It can stabilize the electroweak scale against the dangerous radiative 
corrections. Here we briefly discuss the SUSY extension of the family gauge replicated 
model and also its non-thermal leptogenesis. 

In the family gauge replicated model with SUSY, we have to introduce two Weinberg- 
Salam Higgs (chiral super) fields, say 0„ws and 0dws; in order give masses to both the up- 
type quarks and the down-type quarks, respectively. Further, to ensure the cancellation 
of gauge and mixed anomalies, we introduce mirror (chiral super) fields of our Higgses, 
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say, W, T', p', u', and (p'^g, for which the gauge charges are as if they were complex 
conjugates of W, T, p, uj and 0ss- Then, in the SUSY models, the number of the (chiral 
super) fields for Higgs becomes doubled. 

The first question we would like to ask is: "Can we reproduce the successful fitting of 
the masses and mixing for fermions which we performed in section 3?" . The answer is yes, 
to the first approximation. From the gauge invariance (z.e., six f/(l)'s) we can write down 
the fermion mass matrices as Eqs. (0), (p3D, (|I6|) and (0) by replacing (0ws) {(pdws) 
and ^0ws) ~^ (0Mws)) and also by replacing the dagger Higgs fields to the mirror (chiral 
super) fields, e.g., W"^ W, as we must do since the superpotential shall be holomorphic. 
If we take {4>uws) = {<Pdws), the overall scale of the Dirac matrices is the same for the 
up-quarks, down-quarks, charged-leptons, and Dirac neutrinos, but is suppressed by l/\/2 
comparing with the non-SUSY model which is outside of our discussion. Furthermore, 
from the D-teim flatness conditions, we expect that the VEVs of our Higgs fields are the 
same as those of mirror fields, e.g., (W) = (W). It should thus be noted that in this 
use we have only the five VEVs (W) = {W), (T) = {T'), (p) = {p'), {u) = {uj') and 
(0ss) = (0ss) the fitting even though we had introduced the doubled Higgs fields in 
the model compared with the non-SUSY model. Therefore, the structure of the Dirac 
and Majorana mass matrices remains unchanged at the cutoff scale, taken as the Planck 
scale. If we neglect the effects of the evolution by the renormalization group equation, we 
obtain the similar results given in Table 2. Therefore, to the first approximation, we could 
say that we keep the successful fitting of masses and mixings even with SUSY introduced 
into the model. 

The non-thermal leptogenesis is found to work in the SUSY version. By introducing 
SUSY, the total decay rate of the inflaton through the interaction Eq. (0), the decay 
rate asymmetry parameters of Majorana neutrinos e^, and the degree of freedom of the 
thermal bath of the universe change only within a factor two which is also out of our 
discussion. Then, the inflaton of mass ~ 2M2,3 induce the reheating temperature 
of Tr ~ lO^GeV and hence the correct order of baryon asymmetry is produced via the 
proposed scheme. Interestingly, with such low reheating temperatures we can avoid the 
cosmological gravitino problem [p6| . 

Finally, we would like to ask; "Can we suppress the interactions between the Weinberg- 
Salam Higgs fields and the inflaton?". The potentially dangerous interactions would be 
induced from the term in the superpotential as W = (t><t>uws<t'dws- This term can be killed 
by using the chiral U{1)r symmetry and the global U{l)x symmetry. We assign the 
i?-charges for the Higgs fields and the inflaton, and 1 for the matter fields. On the 
other hand, we assign the X-charges for our Higgs fields {W, T, p and uj) and the 
inflaton, — 2x for two Weinberg-Salam Higgs fields and the seesaw Higgs, and x for the 
quark and lepton fields (x is the charge corresponding to U{l)x)- This symmetry forbids 
uws'Pdws as well as the /i-term W — p(f)uws4'dws- The required /i-term may be 
generated by the effect of the SUSY breaking. Therefore, we can avoid the fast decay 
of the inflaton into the Weinberg-Higgs fields. Notice that the dangerous proton decay 
induced by the dimension five operator is also absent due to this f/(l)_R symmetry. 
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As we had briefly sketched, in the family replicated gauge group with SUSY, we have, 
to the first approximation, the successful fitting of masses and mixings for fcrmions as 
well as the non-thermal leptogenesis via heavier Majorana neutrinos produced by the late 
decays of the infiatons. 

5.2 Further generalization, abstracting the model 

It should be remarked that the idea of how to get the heavier seesaw neutrinos produced 
non-thermally in decay of an inflation boson, could be generalized to letting the boson - 
replacing the boson above - be any sufficiently long lived boson at the end decaying 
dominantly into the (heavier) seesaw neutrinos. Then it is possible that at the stage of the 
cosmological development when the temperature passes the mass of these particles, 
they become non-relativistic and we obtain a matter dominated era. Such an era will end 
when Hubble expansion rate, H, becomes of the order of the width by the decay of 
these particles - pretty analogous to the "reheating" talked about above. 

That it is a true inflaton does not really matter, we could get by the same calculation 
as wc did the same leptogenesis anyway. In this way we could escape completely mass 
estimates suggested for the inflaton. However, since the inflaton mass is very strongly 
dependent on the form of the inflation effective potential used, there is really no general 
mass estimation for the inflaton (in a model independent way). Therefore, it was correct 
that we used the mass for fitting m<^ ~ 2M2,3 ~ 10^° GeV. 

6 Conclusions 

In this article the leptogenesis by the decays of the heavier Majorana neutrinos N2 and 
A^3 was considered. We pointed out that the N2 and A^3 decays could be responsible 

for generating the baryon asymmetry in the present Universe, if they are produced non- 
thermally in the inflaton decays. In our picture then the reheating temperature of the 
inflation is sufficiently low that the wash-out processes mediated by the lighter Majorana 
neutrino are ineffective. 

Accepting this picture for baryogenesis would release the family repUcated gauge group 
model from its only severe deviation from agreement with experiment - the amount of 
baryon minus lepton number produced in the seesaw era - and thus make it agree within 
its pretended only order of magnitude accuracy with all quantities it predicts. 

This means that with the combined picture of an inflaton as suggested here and the 
family replicated gauge group model - in the latest version - we have a very viable model 
where we used the five parameters to produce the fit of Table 2 and some adjustment of the 
inflation properties essentially all the parameters useful for going beyond the Standard 
Model, however, only order of magnitudewise. Indeed we have now fit all the quark 
and lepton masses and mixings - including the neutrinos and the CP violation - and 
the baryon asymmetry, and further successfully coped with bounds on proton decay and 
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lepton flavor violating decays. Also bounds on the neutrinoless beta decay are respected 
though so that neither the positive findings of the Heidelberg-Moscow collaboration of 
(I'^eel) ~ 0.39 eV nor the LSND neutrino are welcome in our scheme. 
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